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Abstract 


This paper demonstrates the invariant distribution of a stochastic dynamical system. We 
give the invariant distribution and numerical examples. We also present a further discussion 


on the computation details. 
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1 The stochastic dynamic system and invariant distribu- 
tion 


The stochastic dynamic system we focus on is 


dx = vdt (L1) 
dv = —VVdt — «(v — Ax)dt + Avdt + edo, 


where i 
V = 57 Bm. 


This equation can be written as 


where 


0 I 
M= f 
er nn 
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We can obtain the solution of the initial form of the problem is 


i 0 0 
a(t) = eM'a(0) +f ett) dwz. 
0 0 I 


Remark. We assume that 


lim eM* — Q, 
i—-4-oo 


Then z(t) will converge to a limit distribution when t > oo. It can be proved that it is a 
normal distribution M (0, C), where 


ge 0 0 T 
C=0 exp(Mt) - exp (Mt) dt. 
0 


'Therefore, we obtain a linear system below 


uc cut - e f 
0 


=o? l| x (eam | : j ) exp am) dt 
= foo 
COT gr) 


Then we can obtain the theorem below 


i (nessun | : : exp (Mt) + exp(Mt) | : 


: exp (M?t) un) dt 


Theorem 1.1. Suppose that we have an equation 


dx = vdt (1.2) 
dv = —VVdt — «(v — Ax)dt + Avdt + odw. 


0 I 
yA— B A-yI 


Mt L0 


Let 


E 
Í 


Suppose that 


lim e 
t—+00 


then the limit distribution is N(0,C) where 


MC+CM!? = —g? of 
0 I 


, 


From this equation we can not only solve the matrix C easily but also obtain some corollaries. 


Here we give the most obvious one. 
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Corollary 1.1. Suppose that 


Q 
II 


Cy C2 
C3 C4 ]. 
We have 

(1) C2 is a skew-symmetric matriz. 

(2) the diagonal elements of C» are all 0. 


Proof. From 
0 


MC --CMT = -0° E 


We can obtain 
C» + C3 = 0. 


But C is the covariance matrix of a normal distribution, so C is symmetric. Therefore, 


C3 = CT 


We conclude 
So (1) is proved. (2) is a corollary of (1). 


2 Numerical examples 


We give two examples. 


Example 2.1. 
0 b 
n=2,0=1,y=1,A 
0 0 
0 0 1 0 
M= 0 0 0 1 
-1 by —y b 
0 -1 0 -w 
C11 C12 0 C14 
cul %2 e -e 0 
0  —ci4 €33 C34 
C14 0 C34 C44 
0.0 0 0 
D= 0.0 0 0 
0 0 -1 0 
00 0 -1 
MC+CMT=D 


We obtain 


1 b 1 b 
ai=] (30? -- 2) , c12 9:02 » C33 (b? + 1) , c34 C44 = 5, C14 i 
Therefore 
+ (3+2) 8 0 —$ 
b 1 b 0 
C= 2 2 4 
0 à 3(% +1) 3 
b b 
zu 0 4 3 
p= zae- 2 (q1, q2, P1, P2) c! (q1, 92; P1, P2) r) : 
From the Fokker-Planck equation [1] 
u o OV g^ p 
e (bip) Opi ðq; 2a A + » Aijpj p st 2 ap op ^ 
we can confirm this result. 
Example 2.2. 
n=1A=a,B=b 
0 1 
M lI 
ya—b a-y 
From 
0 0 
MC+CM? = 
0 —o? 
we can obtain " 
C= | Mayra 
~ 2(a—-3) 
—_ 1 1 -1 T 
P= asc =P (— 36a sc (o9). 
From the Fokker-Planck equation 
L o OV g^ o? 
ax (pip) Op, ðq; (C2 Av + » Aip; | e| +> 2 8p op” 


3 The computation analysis 


In this section, we take more discussion on 


MC - CMT = -0° a 
0 I 


202304.01048v1 


chinaXiv 


Assume that C’ = iC, and in this case, we obtain: 
0 0 
MC' c C'MT « — 
0 I 


For conveniently writing, we denote C" as C in the following 


MC --OMT =- me ' 
0 I 


where M is known, C is unknown based on simple calculation, we have the conclusion that the 
total account of the unknown variable is 2n x 2n. The system is equivalent to a linear equation 
set DX = y, where d € RIXA? y e g^, y € qj". Since 


C= Ci Cy . 
—C9 C4 
where C is a symmetric matrix, C2 is a skew symmetric matrix, C4 is a symmetric matrix, the 


amount of the unknown components, in fact, is only 


Then we obtain that 


MC--CMT = 


7 -C3 C4 [P+ OR 

~\ PCL-QC PC+QC4 Ca —C2PT AGO" 

0 C4 + Cip? + CQT 
C4 + PCi -QC BO TOU -OUPT + CQT 


where P = yA- B,Q=A-qI. 
Thus we obtain that 


Ca + C,PT + CoQ? =0, 
| 4 1 2Q (31) 


PC + QC4 — CoP? 4-040" = T. 
According to the first line of (3.1), it holds that 


D = -C PT 00. 
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Besides, we know C4 is symmetric matrix, which means CT = C4 Thus we obtain: 


Q, =—C,PT — CoQ? 
=- (a eo 
= -PC - QC? 
= —PC, - QCr. 


(3.2) 


Substitute (3.2) into the left hand side of the second line of (3.1), and then we obtain that 


PC2 +Q (—C,PT) — CPT +(PC1 + QC2) QT — —I, 


and 
PO, — POQT — CaP? —QC,PT = —I. 


Thus calculating the covariance matrix C which is equivalent to solving (3.1) is finally equivalent 


to solving the equations 


Ca + CPT +Q =0 
(3.3) 


PG 09^] FAO-a PY — 7 
without the integral. 


Example 3.1. 


so 
—-] b —-1 b 
p- i 
We denote C, by Ch, C2 by C7; and C4 by C;;. According to symmetry or skew-symmetry of 


ij^ 
C4, C5, C4, we obtain the linear equations: 


2C1, + 20€]; — 2C? + 20 (Ch — bCh + C2) = -1, 
Cl; Ch | bC 29 Ci Ch Fb (Coe | C32) =0, 
Ch bC is Ch, bCi, Ch = 0, 

Cá bC} Cå bC Ch = 0, 

-Ci = Ch E Ch =0, 

-C3 z C2, g Cy =0, 

—20L, - 262, 2 —1, 

ch = Cá — 0, 

Ci; = Ch = 0, 


202304.01048v1 


chinaXiv 


Ch + C3, = 0, 


C = 0, 
C2, = 0. 
Solving the equations, we obtain C as 
e+) $ o = 
b 1 b 0 
C= 2 2 4 
0 zg alir) | 
b b 
-1 0 4 3 
Fokker-Planck Equation helps us to check the solution 
dp 7 Ó o OV o? o? 
= i i Aijqj Aijpj — Awe = 0. 
2 2. ag, PP) 7 Bp, a P 2 jd +h spi | e| +> ag 


4 Conclusions 


This paper demonstrates the invariant distribution of a stochastic dynamical system. We give 
the invariant distribution and numerical examples. We also give the details of the computation 


analysis. 
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